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abstract 


The cumulative damage of aluminium alloys with different 
yield strength and various ductility due to seismic loads has been 
studies. The responses of an idealized beam with a centered 
mass at one end and fixed at the other end to El Centro's and 
Taft's earthquakes are computed by assuming that the alloys are 
perfectly elastoplastic materials and by using numerical tech- 
nique. Consequently, the corresponding residual plastic strain 
can be obtained from the stress-stra in relationship. The re- 
vised Palmgren-Miner cumulative damage theorem is utilized to 
calculate the fatigue damage. The numerical results show that 
in certain cases, the high ductility materials are more resistant 
to seismic loads than the high yield strength materials. 

The results also show that if a structure collapse during 
the earthquake, the collapse always occurs in the very early 
stage. 
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INTRODUCTION 

There are two different approaches to structural design. 

One is based on the concept of allowable stress and elastic 
behavior of materials and the other, on the concept of ultimate 
load and inelastic behavior. Since the structure designed by 
using plastic theory is more economic than those by applying 
elastic theory, plastic design concept has been gradually accepted 
by most engineers and has been widely applied in constructions. 

If the external forces are large enough, the plastic- 
design oriented structures, undoubtedly, will undergo plastic de- 
formation and create residual plastic strain. The residual 
plastic strain will gradually accummulate if the structure is 
subjected to dynamic loads such as wind force, seismic loads, 
etc. If the cumulative plastic strain reaches the critical 
values* the structures wi 1 1. col laspe. This process is usually 
termed low cycle fatigue failure. Hence the study of low cycle fatigue 
failure is of importance to the structures designed by using 
the concept of plastic theory. 

There are many papers dealing with inelastic behavior of 
materials (2, 9» 10, 12, 13? and the fatigue damage of structures 
under the action of seismic loads. (4, 5, 6, 10). Few of them 
have stressed the importance of the comparison of different fatigue 
life among the same materials with different yield strength and 
ducti 1 ity. 

*The numerals in parenthese refer to the list of references. 
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Since low cycle fatigue damage can'be expressed in terms of 
plastic strain, the fatigue life of structural systems can always 
be improved by choosing the materials with high ductility. However, 
high duct i 1 e materials general 1 y have low yield strength wh i ch 
will cause the larger deformation than the materials with high 
yield strength. It becomes a compl icated problem to choose 
material for structural members to resist dynamic loading. Based 
on the given* random behavior of plastic deformation, a methrod 
for evaluation of alternative materials to be used in structures 
subjected to random White Noise excitation was presented (15). 

Since there is not any white noise excitation in reality, it is 
desirable to study the fatigue damage of materials subjected to 
earthquake laod in more details. 

The purpose of this investigation is to study the fatigue 
damage of materials with different yeild strength and various 
ductility due to seismic loads. The revised Palmgren-Mi ner 
cumulative damage theorem expressed in terms of plastic strain 
is used to compute the fatigue damage. As an example, aluminium 
alloys (5052 .and 3030 groups) are considered on an ideal beam, 
subjected to 19^8 El Centro's and 1951 Taft's earthquakes. 
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LOW CYCLE FATIGUE DAMAGE MODEL 

In computing the fatigue damage, the Palmgren-Miner cri- 
terion (8), because of its simplicity, has been most widely used 
during the past decades; This theorem can be mathematically 
expressed as follows: 

N 

X. d; = D, <1> 

i = l 

Where d; is the fatigue damage during its cycle, and D is the 
total cumulative damage under N cycles. The theorem states 
that the fatigue failure will occur if the total damage D reaches 
some critical level, which depends on material properties. 

At present, this study is concentrated on the low cycle 
fatigue life which usually is referred to the case with life 
less than 10 cycles. If a structure is supposed to collaspe 
within such a limited cycle, the magnitude of cyclic loading may 
be large enough to cause the plastic strain in the materials. 
Hence the Palmgren-Miner criterion for low cycle fatigue life 
was developed to associate with plastic strain. Manson 
and Gross et. al. (3, 7) proposed the following low cycle fatigue 
life criterion by introducing plastic strain for predicting the 
total number of cycles of reversed-stra i n to cause fatigue fail- 
ure: 

N m (Ae t ) = C <2> 

Where £ e t is the plastic tensil strain, m is a constant depend- 
ing on material properties. N is the fatigue life and C is 
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some constant. 


Later Yao and Munse (14) suggested that, for 
uniaxially loaded metal, the cumulative damage can be expressed 
in terms of material ductility and plastic strains as follows: 


N , ( A e-t- ) 1 /m y — 1 

1 e f 

i ■ 1 


< 3 > 


Where ef is the ultimate plastic strain and is constant for a 
given material, Ae t is the plastic strain at tensil cycle and 
m is material constant dependent on ambient and loading rate. 
Usually m is defined as function of plastic strain in tensile 
and compressive cycles, denoted by e c and et, respectively. 

1/m = 1 - 0.086 ( A et-) < 4> ^ 

A e c 

The cyclic history of plastic strains in tensile and compressive 
cycles is illustrated in. Figure 1. Equation <3> and <4 > will be 
considered as the fatigue damage criterion through this invest- 
igation. 



Figure 1. Cyclic History of Plastic Strain 
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MECHANICAL MODEL 

Consider a beam with a centered mass at one end and fixed 
at the other end subjected to a fluctuating excitation as 
shown in Figure 2. 



Figure 2. Mechanical Model 

The equation of motion of this spring-dash-pot system is 
well known. 

M x + C x + G(x) = M Y <5> 

Where M is the concentrated mass of the beam. 

C is the damping of the mechanical system. 

x is the deformation of the beam. 

y* is the accelerogram of seismic loads. 

G(x) us tge f i rced-deformation function dependent on the 
stress-strain curve. 

The force-deformation, function, G (x) , can be be derived from 
the stress-strain relationship of a given material. For sim- 
plifying the . computat ion , the aluminium alloys considered here 
are assumed to be elasto-plastic materials. The stress — 
strain diagram of such a material is shown in Figure 3. 
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Stress 



Figure 3. Elasto-plastic stress-strain diagram 
Let A be the cross-sectional area of the beam. 

E be the modulus of elasticity. 

L be the length of beam. 

e r be the residual strain. 

e y be the yielding strain. 

fyt be the yielding tensil stress. 

fy C be the compressive yielding stress. 

If f y -t = fyc = fy,' the force-deformation function can be formulated 
as follows: 

G (x) = EA ( X - ej. ) If I X “ e r I < I e y| 


G (X) - Afy 


If 


L 

X - e r 
L 


> | e y 


<6> 


Since the analytical solution of equation (5) is not available, 
Ruge-Kutta Numerical intergretion was used to obtain the defor- 
mation, x. From the known deformation and stress-strain curve, 
the cyclic history of plastic strain can be found, then the 
fatigue damage can be calculated from equation < 3> • 
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NUMERICAL EXAMPLES 


In order to illustrate the low cycle fatigue damage of 
structures subjected to dynamic loads, several aluminium 

alloys are studied. The mechanical properties of these ma- 

• , 

terials are tabulated in Table (1). 


Alloy 

Yield 

strength 

ksi 

Ducti 1 i ty % 

Modulus of 
Elastici ty 

10 3 ksi 

Alclad 3003 - 

H 1 2 

18 

20 

10.0 

Alclad 3003 - 

H 1 4 

21 

16 

10.0 

Alclad 3003 - 

H 1 6 

25 

14 

10.0 

Alclad 3003 - 

HI 8 

27 

10 

10.0 

5052 - 

H32 

28 

:i8 

10.0 

> 5052 - 

H34 

31 

14 

10.0 

5052 - 

H36 

35 

10 

10.0 

5052 - 

H 38 

37 

8 

10.0 


Table 1. Mechanical Properties of Aluminium Alloys 


Since the natural frequency can not be well defined if the 
structural systems allow to undergo plastic deformation, the 
terminology of "artificial frequency" is adopted. Let p b’e 
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the . arti f i ca i 1 frequency which is defined as fl AE) / (ML). If 

both sides of equation <5 > , are divided by the mass M, it becomes 

i 

x + 2 c p a x + g(x) =*y < '7 > 

Where c is the damping ratio, p a is the artificial frequency 

and g(x) is G(x)/M. It is to be noted that g(x) can be reformulated 

in terms of artificial frequency. 

In order to make the comparision, two seismic loads are 
utilized in this study. One is 1 948 El Centro's earthquake 
NS component with maximum peak in the order of 0.312g, the other 
is 1951 Taft's earthquake S69E component with a peak value 
in the order of 0.157g, where g denotes the gravitational force. 
Because of the high peak value of El Centro's earthquake, the artifi- 
cial frequency of the mechanical system varies from 50 rad/sec to 
7Q rad/sec, and for the case of Taft's earthquake, the frequency 
varies from 35 rad/sec to 45 rad/sec. With the length of beO, L, 
being 10 inches and the information shown in Table 1, the fatigue 
damages of aluminium alloys subjected to El Centro's and Taft's 
earthquakes are computed by the aid of 360/65 IBM computer. The 
computer program is listed in Appendix. The results are tabulated 
in Tables 2, 3, and 4, and also plotted in Figures 4, 5, and 6. 
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Table 2 Fatigue Damages of 5052 Group Alloys Subjected to 
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El Centro's Earthquake 



Figure h. Fatigue Damages of 5052 group 
alloys due to El Centro's Earthquake 
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Figure 6. Fatigue Damage of 5052 Group Alloys due to 
Taff's Earthquake. 



D I SCUSS I ON AND CONCLUSION 


For Table 2, 3, and k or Figures 4, 5, 6, they show that 
in general, the fatigue damage will increase as the stiffness 
of the beam decreases. However, in some cases, the fatigue 
damages of higher stiffness member are smaller than that of 
lower stiffness member as ic can be seen in frequencies 62.5 
rad/sec and 65.0 rad/sec for a] loys .5052-H32 , and 5052-H36. 

This phenomenon is dur to the fact that the frequency of the 
beam coincide with the frequency of seismic load. 

Since a structural system is usually defined as "collapse" 
if the cumulative fatigue damage reaches one, the discussion 
here will concentrate on those. materials whose cumulative 
damage smaller than one. Under Taft's earthquake, the cumula- 
tive damages increase as the yield strength of alloys 5052 group 
decreases regardless the ductility of materials. But it is not 
the case for the same materials subjected to El Centro's earth- 
quake. The results show that in certain frequency domain, the 
material with high ductility is stronger to resist seismic load 
than the material with high yield strength. For easy comparision, 
the cumulative damage historigrams are plotted versus time as 
shown in Figure 7 and 8 for alloys 5052 group in the frequencies 
of 60.0 rad/sec and 8 for alloys 5052 group in the frequencies of 
60.0 rad/sec and 62.5 rad/sec. 

The damage historigrams are also plotted for alloys 5052 subjected 


15 = 




to Taft's earthquake in some frequency domain. It is very inter- 
esting to note that the cumulative damages pile up very fast 
in the early stage of earthquake action as shown in Figures 
7, 8, 9, and 10. In other words, fatigue damage is always 
created during the first 5 - second period for the 30-second 
duration El Centro's earthquake, and the first 12-second period 
for 60-second duration Taft's earthquake. It can be concluded 
that if a structure can survive during the first quarter period 
of earthquake action, it surely will not collaspe and under this 
earthquake. On the contrary, it can be said that if a structure 
collaspses during the earthquake, the failure always occurs in the 
very early stage. 
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Figure 7. Fatigue Damage Historigram for 5052 Group Alloy's in the Frequency of 60 rad/sec 
Due t*b El Centro's Earthquake. 
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gure 10. Fatigue Damage Historigram For 5052 Group Alloys in the Frequency of 
42.5 rad/sec due to Taft's Earthquake. 
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10 = 01 ^ :• ■ ' '! ’ 

1 io*oi>o. 5 ■ ;v. ... . f .. : _ 

IF (ID-IID) 80,90,100 ' . :. f ' \ i 

_ I D= IJCLt.l — ' - .. . - 

GO TO 90 . •- 

WRITE C 6,904) I 

FORMAT f 2X, 12HERR0R I N AC ,13) ‘ ' . 

id=io-i 1 , ; ..... . j 

IR = IRMD ' ... • X" 

_ACC£L_U 1^1).=ACIJJ.»0«386.4 : ; 

00 120 J=J1, IR i v 

ACCFL ( J ) = AC C EL < J- 1 ) .♦. C AC 1 1 *1 IrAC.Ci I i *0. 386.4/0 L ..i 

J 1 = I D+J 1 

CONTINUE ■ ' ■ _ 

READ! 5,907) APH,BPM,CPM 

FORMAT! 3F 10»6|* L ’ 


READ! 5 , 90C ) N 

_MDIM=N*2 ; _ 

FORMAT! I 5) 

WR I TE ( 6, 999 ) ' ; 

FORMAT Cl HI ) 


- — •- — ■ ■ , J ••• 


• jl ; 

A 1} A ( 


• ■ < • • •- 
; , i.iL 5 ..a.u. 



LEVEL ?'> 


i 


MAIN 


DATE » 72215 


17/18/44 


„3i READ l 5,932) ...( TL(.I 1,1 =UNL 

90 2 FORMAT! 2F10. 51 

PF\0(5,o08) TLMU, LCHECK 

908 F 0 R M A T ( F 10.5,15) 

DO 20 1*1 f N ..... . 

•'RCAn(5.O0H OMGAm.ETAm .DULTVm »EYELDm *CON(n 

901 FORMAT! 5F10.5) ; 

YCLD ( I )=EYFLD( I )*Tl< I) 

0! 11=0. J 

YDMIN! I 1=0. 

Ycmxm=o. 

QMF ( I 1 =DMGA (I 1**2* YE LOCI) 

SKI! I 1 =OMGA ( I 1**2. ■' ‘ 

S K 2 ( I 1 =SK1 ( I lAf.ONlI 1 

OMAX ( T 1 =QMF ( I 1 ' 

omin! i >=-ome( 1 1 •*. 

YMAX! I 1 =YELD( I 1 ' . * ‘fj /'*N 

YMIN! ! )=-YELD{ I I ■? - 

WRITE (A, 923 1 OMGA ( I 1 ,ETA(I 1 .OULTY! 11 .YELP (I 1 ,CON! I 1 »TL! I ) 

923 FORMAT ( 2X , 10HFREQUENCY F10.5.17H DAMPING RATIO Flb.5»ilH DUCTILI 



NL=NT 

CALL DAMAGE! YD I SI, YE LD( 1) yDAMGCJ 11, 6 AMG V ! JJI , l< { 1 ) , KT I MJL, TL ( 1 ) 9 

1DUL TYl 1 1 , NL ,KT1 , 1 1, ~ . - - -• 

CALL DAMAGE! YD I S2 , YE 10(2 ! yDAMGC ! 2 1 , DAM GV f 2 ) , K! 21 *KTIM2,TL« 2 
1DULTY! 2 1 »NL ,KT2 , 2 1 ‘ 1 ' , a.. ... 

TFND1 = XT(KT>M11 " ** ‘V'”’ 


TEND2 = XT!KTIM21 * 

WRITE! 6, 800 1 TFN01 ,TEND2 . _ _ 

800 FORMAT ( 2 X , 3HT 1 = l P E 1 6 • 7 » 1 b X f 3HT 2 = 1 PE 16. 7 j 
DO 30 1=1, N 

WRITE! A, 9211 N, YDMAX! I ) , YDMIN! J 1 
30 WRITE 16,920) N , DAMGC 1 1 1 , DAMGV! I) , K Uf 

920 FORMAT! 1 OX, 13, “6H STORY, 16H DAMAGE FACTOR 1P2E16.7»20H NUMBER OF 

1 CYCLES 17// 1 \ 

921 FORMAT !iOX, I 3»6H ST0RY»19H MA XD ISPLACE ME NT 1 PE 1 6« 7 9 19M ~ MIN DISP 
1 LACEMENT 1PE16.7//5 


K1 = K! 1 





nn i'»o ! = l, ki 

J =KT l ( 1 ) 

160 WIUTF(6,9A0) I,XT(JI 
no IS* I = i , K2 
J =K T 2 ( I ) 

150 WIUTEC6.950) .I.XT1JL — .... 

9 40 FORMAT ( 2X ,1 2HNO OF CYCLES 15 ,2X , 5HT I ME = 1 PE16. 7,2X , 2H0 1 I 
950 FORMAT ( 2X» 12HNO OF CYCLES 15 ,2X * 5HT1 HE* 1 PE16 . 7 ,2X , 2H02 I 
IF(LCHECK), 31»32»31 

- ...32 CONTINUE 

STOP 

end ; : 


c 

c 


S'lRROlITTNP RKC,S{ P«MT,Y t r)FRY t Nn.TM» IHLF, FCT«nUTP« AIJX » 


nr MF NS TON Y ( ? ) ,f)E n Y ( 2 I t MJX ( 8 » 2 ) ,.A I 4 » , B ( 4 ) * C 14 1, l»RMT t 5 > 

no j i = i. , no I m 

1 A')X(R. T )=.06666 6A7»nFRVm ; 

y = f> n m f ( i ) 


XFN r »= n RMT { 2) 

M=PRMT(31 
poMT( S ) =0, 

CALL rCT(X,Y,OFRY! 


C .-F^OHR TPST / 

TF(M*(XFNn-X) ).3R, ’.7,? _..J 

0_ PRF°APATTPNS FOR RtJNGE-KU TTA MFTHOrV* U‘ . . 


. .. 

n - , 


2 At n = .s 

AL?.>_= . ZJ1Z3J112. 


A ( X )*1 .707107 
A { 4 1 =. 1 666667 
Mil:’, 

.*.( 2 1=1. 

_?.L4.L=2-* 


-r 1 — : 




V 


.V 


TV* 


»•“/ . »i', ! 

>*■- — -r* P — ■ : 


-Oi 4 . 

A 

i ( A 


* • \ 


ULL 


t 


r m = . s 

r ( = 

C ( ^1=1.707107 
r ( 4) = . 9 

PREPARATIO NS OF F 


v - ;j ... ,. : ' 

r ; - ^ j| jV.JJ 5 A r »; 


M 
.1 - 


nn 3 l = i 9 NDT M 
AIIX( 1 , I )=YU ) 

A')Y(? , I 1=0FRY( T I 

a 1 1 X ( T 1 = 0 * 

3 AtJX{6»n=0. 

i p rr=o 




H=HfH 
JML c r“l._ 
I STCp=fl 
I r ; NO = n 


S 4 .- : 




START OF A RIINOF-KUTTA STFP 
6 T F (( X+H-XFNQ )»H)7«6«5 

5 H= XFND-X ' 

6 T FN0»1 . 


i 

! . . 


9A, 



FVH. 


720*0 


1 G/22/T5 


' » • ’ » . ' ' ' i ' r i | i f i t r\ t | 

f c ( OR MT | H | ) ^ O t R f /w O 

q I Tr<;T=0 
9 j STFP= I 5T P P+1 


r ART OF INNERMOST RIJNGF-KUTT A LOOP 

n j. 

J = MJ) 

I = P ( J ) !■ _ 

KUI 

) n. l = l.M niM 

= H*DFRY( T ) ' t 

’ = A J ♦ { R 1 -R J* AUX (jh t I).) 

n=Y< i i +Rj> ” ; 

' = R2 + R 2 + R 2 __ o^Tl/jL 

IX ( 6 , I ) = Aijx ( f» , i ) +R 2-C J*Ri Vk v * 

(J-M 12, 15,15 




^ % 




H=.5*H - '-» if • v 

nn 19 1 = 1 , NO 1 M •••'••■ :■;■■■ iV ’ 

• P K 


D K 


Y ( T ) = AIIX ( l , I 1 

RK 


D p R Y { I ) = A UX ( 2 » I) 

RK 

1 9 

AI.JX(6, T ) =AUX(3^n \ 

RK 


GOTO 9 

RK 


• \ '/ 

RK 

r. • 

IN GASP I.TFST=1 TESTING OF ACCURACY IS POSSIBLE 

RK 

20 

[ MO0= I STF.P/2 

RK 

• 

? F ( I STEP— I MOD— I MOD) 21*23*21 f 

RK 

21 

CALL FCT ( X * YfDERY ) 

' t‘ ’ 

RK 


.27 




G LFV C L •>!) 


OATF a 72080 


16/22/35 


DC ..22....1.r 1 • MD.I !1 

A U X ( 5 » I » = V C T ) 

?2 A ijx ( 7 * n=Df?P.Y{ T 1 

r-ITfi q 


fOMOiiT A r f ON np t p S T valuf delt 

Z-UD.C.LXko.i ; 

nn = 1 , MO I m 

.26 DFLT = DCLTf AiJX(P , I ) * A8S ( AUX lA. I i-YII ) )_ 

I c IDFL t -P»MT (4) 1 28, 28,25 


FPPOn IS Ton GREAT 


26 DC 27 1 = 1 .NDIM . ' , 

71 AlIX » I- ) sAUX ( 5 i T 1 _• 
I <;TFP = | S t FP*TSTFP-4 

X = X-H ,_1. 

SfFMO = 0 ■ 


“:;•»***•**■*»**',»* 

.i 1 '* _ •»* •. 


...TV* — 


RESULT VALUES ARF (JOQO s/.. 

2 8 CALL FfT{ x, Y,0ERY) # [ft 

on 7 Q f = 1 v NO T M 4 ? ->••// gf/ jS 

mixu.t )VyuT" " ’ 

ai»x ( ? . T l = r>Fp y t n ; f '/ f 1 

Auxn, n=Aux( 6 ,n w jLC. 

Y { I.) =AUX( 5 t T I. ,.„_2 vJi ±Uii:a :,V.i !lV~ 

2D nFRY( I )=AUX(7 t n • •- 4 - 

c ALL nt.lTP < X-H , Y » DFP.Yj 1 Hl.^.0 I!!kP8MTJ._„;/_ 
IFtPPMT|S|T*Tl,30,40 v\f:-V v,' ' ; 

™ DO 71 Irl.NOI M V. '• /:•; 1 -1 

v i i i =aux ( l ♦ T ) 'j * ‘ A 

?i._ ocry m =aux(. 2.,.U y : - 4 .' . 

TPPC = IHLF ■ % 


Co 


y 

•r . , 


>-.A I ' A b 

— ;V:.l lv 


V / 


. ..r V 


’ \ ii — £ • •• r. 

\ ■ >.*. ,, ; v < R 

I NCR F HFNT GETS DQURLFD \ . A- ^ R 

3? T HLF = { HL F- 1 
TS t FP=TSTFP/2 

- f »•, V- , t'i'. ‘ D 

' ...,nv« R 

H = H + H 

TF( IHI F 14,33.3? 

• ' R 

R! 

3-J I MQn= I STFP/2 

f F { T STFP-l MflD- IMOOl 4,34,4 

\ ' Rl 

^ ?J 


TF(DFLT-.n?*{>RMT{4} 935*35»4 
IHLF = IHLF-l ’ ' 


l STFP=fSTFP/2 
H=H*H 


GOTO 4 


! 

. '• i 

: : — •- -- > 


• N3 
00 

/ ,-v*- 






0 LFVR. . ?n 


PK OS 


HATE o 720 PO 


16/22/3 


_C_ 

r 

c 


*r T iJ*NS TP.. CALLING PROGRAM, , 

^ !Ml.i: = i| 

..TALL FfT ( X ,.Y. f OFR.Y ) ; ' ..... 

GOTn 3 Q 

37 THLF = I2 .... 


r.nTn iq 

..30. I HL P .* 1 3 ; 

3 q CALL OUTP< X,Y,6FRY, IHLP«N0tM» PRMT) 
AO P p T.U R N i 

FND . 


;..f 









U- 


SURRCUJTINE FCT ( X » Y »DEKY } 

DIMENSION Y ( A ), OERYIA )* ACCEL ( 3000 ) * 
10MGA(?) t CTA(2),YMAX(2),YMlN<2),Y2LD(2)»OMAX(?J »QHIN( 2 ) «Q I 2)' »SK1 ( 2) 

2 » SK ? ( 2 I » QMF ( 2 ) . 

COMMON/ VP/YM AX » YM I N » YELD f QMAX » QMI NvQ * SKI »SK2»QM£ 

CCmmoM/.V1/ACC£L...DMGA».EIA 

XN = X/o. 014-1. 

N X = X N '. , . _ ........ 

T F MP =NX 


CALL PI.ASI QFltYI 1) » Y M A X ( 1.1 , YMIN ( 1 ) , YEL0( 1 .» » QMAX U> » QM INUI , Of 1 J « .... 
I5K 1 ( 1 ) » SK2C 1 ) tOME (1) ) 

_ .CAL I PJLA. Si Q F 2 * Y l 2X « Y MAX( 2 1 . YM I N ( 2 ) . YELP C 2 Lt QM AX £ 2 I.QH INI 2 1 a Q1 2 1 a 

1SKU 21 ,SK2(2) ,0MF(.2J I 
DER Y ( 1 ) = Yl 2 ) - . _ ’ 

DERYI 2I=-?.*ETAU)*0MGAI1 1 *Y l 2 > .-OF 1- ACCEL 5 NX ) -*( ACCEU NX*1 i - 
1 ACCEL I NX H* ( XN-XE .MPJ_t? .*JE T A 1 2 5*0 MG.AJ.2.L*.Y.^AJ, + Q ^ 

DER Y ( 3 J = Y t A 5 ■ .v> v ji- . " " ;• ‘-V . 

DE RYI 6)=-A.»ETAC2? »0MGAC2)»YI4>~2o»0F2»2>ETAm»0HG A a ?»Y I 2).»QF1 


RFTURN 

END 


v 


V' !• -'SO:- ;.<- r 

£ X^XX 

( n / 




'• • . r 

...r ;. , /.■ •• 

v / 


A 


»* 

■»:» 



LEVEL 20 


MAIN 


OATE - 72215 


14/09/07 


SUBROUTINE OAMAGE I Y« YELD » DAMGC »D AMGV ,K» KTI ME » TL, DULTY ,NL ,KT , I DX ) 
DIMENSION YOOOO) ,OPP ( 3000) ,DPN< 3000 ), KT ( 3000 S 
K T I M E = 0 
DO 401 I »1, NL 

dpp ( n=o. 

— 401 - DPNI I ) =0 • ; U 

SK=I «/TL 

YMIN=0. 

fpp=o. . .• - : - 

t ■ E PN = 0« •' • 


DAMGV=0 a 
DAMGC=0. 
N = 1 
J = 1 




'zas.:* 


-»■ . V •/ ::«•/ 


!\ ■ ■■ - 


j-i •: “ 

IF (ABSI Y(NI J-YELD) 403,403,404 ^ -Tvi ■ I\r ** ^ . .... 

n=n^i . • j- 

IF- (N-NL ) 402,402V»2 0 , r— , , 4-^ 

IF (YIN)) 419*405*405 J ‘ fu'-*'' .<'* V " - ;V ' 

YYE=Y( N J-YELD — Cl. — _. ' L l’... . 

EPP=YYE*SK * "* v <V> S'l !;- /A Vv;,; 

YMAX=Y ( N I : 1/....^... V-. _ .. .... , 

OPP( J)=EPP-EPN J . //' ■'- ■;. / V \ 

kti j i=n 1 

N=N + 1 h slLv J X • *, / A. . . 

IF (N-NL) 400,400,420 - -- 

IF (Y(N)-YMAXI 410,409,409 u V ' ' t >" > ‘ ‘ 

YYE=YYE + Y(N J-YMAX iXX.. .... » , - 

GO TO 406 . U ‘ a' .v f ■ ' ■ . 

YYN=Y(N}-IYMAX-2^YEL^S 1 ■■ ■ -~ - i — ■ • ■ . •••- — 

I F I YYN ) 411,407,407 ; k 

epn=yyn*sk • -i-' — . - 

YM I N=Y ( N ) • \ v -'. 

I F ( DPP( J ) 1 412,414,413 — - - v*. - - 

WRITE! 6,9991 \ . ' .. • iv ' 

-FORMAT (IOXtI 

RETURN • <*'*’ . 

j _j + j ^ * ' * « ■! ;> ■» ; > ‘ ' 

OPNI J)=EPN-EPP 

N=N+ I .... 1 : 

IF (N-NL) 416,416,420 \ 

- 1 F ( Y I N ) - YM IN I--41 7 ^-41 7 -,4 1-6 X, 

YYN=YYN+Y(N)-YMIN \ 

EPN=YYN*S.K : — * 

YM!N=Y!N! = \ 

GO TO 414 1 — : 1 

YYE*Y8 NS-VMSN-2«4YEL0 


— v - 


•TT-T" 



14/09/07 


OVEL 20. ■' ' DAMAGE . DATE » 72215 


IF (YYE) 415,415,406 - 

419 YYN =Y( N ) ♦YOLD 
GO TO 411 

420 K = J 

00 422 J = 1 » K 
IF (OPPIJH 421,422,421 

421- C0ES = 1.-0.8 6*DPN( Jl/OPPI Jl — 

OAMGV=OAMGVMDPP ( J I /DULTY)**CGES 

IF (KTIMF.NE.OI GO TO 423 - 

IF (OAMGV.GE.l. I KT I ME=KT( J) 


423 OAMGC=DAMGC + COPP( J ) /OULT Y) ** 1. 86 

A WRITE(6,998I OAMGC ,0 AMGV, J, I OK 



SUBROUTINE PlASt QM, Y«YMAX> YMIN, YELD» GMAX V QMI N »Q« SK1 V SK2» QME 9 
IF (QJ 212,201,208' 

201 IF I AOSm-YELOl 202 , 202,203 . 

202 0 M=SK 1 *Y •- ; 

GO TO 216 

-203— IF! Y) 204,206,206 — — * — — 

204 QM = QM INM Y-YMI N) *SK2 . 

205 0=-l. ^ 

YM I N=Y , ■’■••• 


QMIN=QM , • -.1 

GO TO 216 ' 

„QM = QMAX-M Y-YMAXl- 
0 = 1 . ' 

Y^A X = Y 

OMAX=QM 

GO TO 216 . 

IF (Y-YMAX1 209,206,206 

-YMA2=YMAX— Y-2,*YELQ. _ 

IF (YMA2) 210,211,211 

QM=QMAX+( Y-YMAXI* SKI 

GO TO 216 

QM=QMAX-2.*QM6-YMA2*SK2- 
GO TO 205 

IF (Y-YMIN1. 204, 204,2-L3- 
YMI2=Y-YMIN-2.*YELD 
IF ( YMI 2 ) 214,214,215— 
QM = QMINMY-YMIN)*SKl' 

GO TO 216 — . 

QM = QMIN«-2.*QME»YMI2*$K2 

GO TO. 207.. — r - 

RETURN 

END _ — — - 


— 

^JL u 

**»' <* - *■-: 


— — V- \‘->v 

: — ->V i;: L-r- — 

v Kh/ *•-/ x ■ • / \ v •••■ 

■*» X~— r '.^ 1 * t M. ,m ■■■«■ " * !■») ••» '♦*- -• I 


»V H /'■' '*! ' \ / 1 - <r 

il: . u.L*. 4 . LI — 

• > 'Tr'- 1 ; 


* * - * i , * v * * 

/ _ • 

' .. V --- 

• A ;/ . • . 

u.w ; . 


• — t > 4 - 


. * . -o V * , /*. 1 i . * . , 

*■> y [ *. (.i v. « v y :**:t ■ * *“'■ • -l* 


LEVEL 20 


MAIN 


DATE = 72215 


14/09/07 


SUBROUTINE OUTPI X , Y, DER Y , I HLF , NO I M, PRMT > 




• • . V •*’ * f. \ •’ v. . '* • * 

’*• v». ‘ ■ ■ : ; ,. f \ .-V ’ 

•** 1 . • ; : •• • tv- 


DIMENSION Y ( 4 ) , OFR Y ( 4 1 1 PRMTI 5 > YOI SI (3000 ) , YDIS2 ( 3000 i , XT < 3000 i 
l,YOMAX(2),Y0MIN(2l 

- — CCMMnN/vn/YOISl, YDIS2,Y0MlN,YDNAX,NT,XT,DMU,T-?-,APM,BPM,CPM 

N = N 0 I M / 2 

on 212 I =1 »N - — 

J = 2* I- I . 

IF (Y(JJ) 211,212*213 . .... 1 

211 IF(YDMIN( I)-Yf J>) 212*212*214 

214— YOMINI I )=Y(JJ : L» 

GO TO 212 f 

213 IF (YDMAX(I)-VCJ)) 215,212*212 ...._ — • 

215 YnMAX(II=YU| 6 1J +,*"***,„ 

212 ... CONTINUE 

IF ( IHLF- 10 1 216,216,217 ^ ‘ ‘ ./ 

.2 17... WRITE (6, 9 10). X, IHLF. - 

910 FORMAT ( 10X, 2HX= ,F 10, 5»10X, 5H1HLF* , I 3) / i<; ^V\ : .w‘V"* ' • 

.. return i:.\. 

216 IF(X-TT) 227,228,226 ;* A* j\ •• / : >'• ' ' 

2 28 ... XT ( NT ) =X ^ a. J h. - V'.'’ 

YDISKNT)asY(I) v.. * . //..•''* ‘ f \ V X ' ' ' 

YDIS2INT )=Y( 3) 1— ' ■ 1*1 1 

NT=NT+l :;^U \ ■■■,•/ X \ / X : X 

TT = TT + DMU ,JL1, L. U _ 

227' IF (NT-3C00) 225,226,226 "f* \ V. ’ ' 

225 CHEK-ABSIYI 1)) Xi , l 1 1— .... 

DO 219 I=2,N0IM v ; ' / 

.... IF { CHEK-ABS ( V C l i XI— 24-9^24-9, Z l 8 ..L ' .. L~LLL — 1 — 

218 CHEK=ABS(Y(in ; W-. 

219 CONTINUE LI i'J. 

220 IFICHEK.EO.O. » RETURN '-.‘V X • ’ V'* 

- IF (CHEK-l. 0E-.03 ) .. 22 1,22 W222~lv • 

221 PRMT ( 4 ) = APM \ • •/, ''X- *'" v v'" ..V- 

RETURN J — —LL-II ' 

222 ICHEK=CHEK . r ' 

IF ( ICHEKI 223,224,223 , : - : a— I lf! * IL ILi-.l.-J-l t - - 

224 CHEK=CHFK*10. '\ 

... GO TO 22 2 . .. \ : - .' 

223 CHEK I = ICHEK \ 

PRMT ( ^ ) a CHEK/CHEKI*BPM 

P FT URN V 

226 . WRITE ( 6,911'). X — — 

911 FORMAT < 10X*8HNT*3000 F10»5J . 

PRMT l 51*1. _L X-l— : 

RETURN 




•'X;^ -.'v- 

i f;'.v -.v, ... fX.: -r 


>.J- V 'U : 




. \. 


